Derivative Formula and Gradient Estimates 
for Gruschin Type Semigroups* 

Feng-Yu Wang 

School of Mathematical Sciences, Beijing Normal University, Beijing 100875, China 

and 

Department of Mathematics, Swansea University, Singleton Park, SA2 8PP, UK 
Email: wangfy@bnu.edu.cn; F.Y.Wang@swansea.ac.uk 

April 23, 2012 



Abstract 

By solving a control problem and using Malliavin calculus, explicit derivative for- 
mula is derived for the semigroup Pt generated by the Gruschin type operator on 



i=i jh=i ^ 



where a G C^(IR™';IR'^ (8) M"^) might be degenerate. In particular, if (t{x) is comparable 
with \x\^Idxd for some / > 1 in the sense of (|1.5p . then for any p > 1 there exists a 
constant Cp > such that 

\VPtf{x,y)\ < ^M^^Ml^Jh^, t>0,fe ^,{W-+'),{x,y) G 

Vt A y/t{\x\^ + ty 

which implies a new Harnack type inequality for the semigroup. A more general model 
is also investigated. 
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1 Introduction 



It is well-known that a liypoelliptic diffusion semigroup on M has a smooth transition density 
w.r.t. the Lebesgue measure (cf. [IE]). An interesting research topic is then to derive explicit 
estimates on the derivatives of the diffusion semigroup. To this end, the derivative formula, 
which is called the Bismut formula or the Bismut-Elworthy-Li formula due to [H [11], has 
become a powerful tool. 

In the elliptic setting, the formula can be explicitly established by using the associated 
Bakry-Emery curvature tensor. But in the degenerate case the curvature is no-longer avail- 
able and the existing formula established using the Malliavin covariance matrix is normally 
less explicit, see e.g. [H Theorem 10] and [21 Theorem 3.2]. To establish explicit derivative 
formulae for hypoelliptic semigroups, one has to build and solve some control problems as- 
sociated to the corresponding stochastic differential equations, see e.g. [HI [191 EHj for the 
study of generalized stochastic Hamiltonian systems, and see [21 Section 6] for some simple 
examples. See also [IT] for the study of hypoelliptic Ornstein-Uhlenbeck semigroups. 

Among Laplacian type hypoelliptic operators without drift term, two typical models 
are the Kohn-Laplacian on Heisenberg groups and the Gruschin operator on M^. In recent 
years, the gradient estimate and applications have been intensively investigated for the heat 
semigroup Pt generated by the Kohn-Laplacian on finite- or infinite-dimensional Heisenberg 
groups, see [H [SI [T21 [S] and the references within. In particular, the gradient inequality 

(1.1) T,{PJ)<CPtT,{f), t>ojeClt>o 

is confirmed in [12] for some constant C > 0, where Fi is the associated square field. This gra- 
dient inequality has important applications, for instance, it implies the heat kernel Poincare 
inequality and thus (cf. [4]), 

(1.2) Tr{PJ)<-^P,f, feClt>0 

for some constant c > 0. 

Accordingly, one may wish to prove (11. ip and (II. 2p also for the semigroup generated by 
the Gruschin operator + x^^dy on M^, where / G N. As pointed out to the author by 
the referee that when / = 1 these can be confirmed by using the known inequalities on the 
Heisenberg group and the submersion ip : {x,y,z) i— )■ {x, z + Indeed, letting Pt and Fi 
be the semigroup and square field associated to the Kohn-Laplacian X'^ + on M.^, where 
X := dx — ^dx, Y := dy + ^d^, we have 

X{fo^) = {dJ)o^, F(/o^) = (x9,/)o^, feC\R'), 

so that 

Fi(Pi/) o ^ = T,{PJ o ^), {PtT.if)) o ^ = P,fi(/ o tjj) 

hold. When / > 2, (11.11) is however not yet available. We also would like to mention that 
for / = 1, the generalized curvature- dimension condition introduced and applied in [3 [SI [Z] 
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holds, so that the corresponding results, in particular the gradient estimates and applications 
derived in [6], are valid. Even when / > 2, although their generalized curvature condition 
is no longer available, a more general version of curvature condition has been confirmed in 
[IE], so that the L^-gradient estimate as in Corollary 11.21 below for p = 2 holds. 

In this paper, we aim to establish the Bismut-type derivative formula and gradient esti- 
mates for the semigroup generated by the following Gruschin-type operators on 



^ i=i jk=i ^ 



where a G C^(M™; M'^®R'^) might be degenerate. In this general case, it seems hard to adopt 
the above mentioned arguments developed for Heisenberg groups and subelliptic operators 
satisfying the generalized curvature. Our study is based on Malliavin calculus. 
Let Fi be the square field associated to L. Then 

(1.3) T,{f){x,y) = \Vf{;y){x)\' + \a{xrVf{x,-){y)\', (x, y) G M"^+^ / G C^(M"^+'^). 

We will use | ■ | and || ■ || to denote the Euclidean norm and the operator norm respectively. 

To construct the associated diffusion process, we consider the stochastic differential equa- 
tion on W^+'^: 



[lA) 



dXt = dBt, 



where {Bt, Bt) is a Brownian motion on M"^"'"'^. It is easy to see that for any initial data the 
equation has a unique solution and the solution is non-explosive. Let E^'^ stands for the 
expectation taken for the solution starting at (x, y) G M™"'"'^. We have 

Pi/(x, y) = E^'yfiXt, Yt), / G (x, y) G M™+^ t > 0. 

To establish explicit derivative formula for Pt, we need the following assumption. 

(A) For any T > and x G M™, Qt ■= Jq <^{x + Bt)cr{x + Bt)*dt is invertible such that 

Obviously, Qt is invertible if so is cr(x) for a.e. x G M™. According to the proof of Corollary 
11.21 below, assumption (A)is ensured by (11. 5p below. 

Theorem 1.1. Assume (A). For any f G Cl(W^^'^) and v = {vi,V2) G M'"+'', 
V,PTf{x,y) = E-'^{/(Xt,Ft)Mt}, {x,y) G T > 
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holds for 

Mt =^^^ - Tr (g^^ ^{(V.,a)a*}(x + B,)dt"j 

where Vi, stands for the directional derivative along v. 

To derive explicit estimates, we assume that a{x) is comparable with \x\^Idxd in the sense 
of below. 

Corollary 1.2. Let I G [l,C)o) and assume that 

(1.5) ||ct(x)|| > a|x|', ||a(x)|| + ||Va(x)|| ■ |x| < b|a;|', x 

holds for some constants a,b > 0. Then for any p > 1 there exists a constant Cp > such 
that for any v = {vi,V2) & M™+'^, 



;i.6) \V.PTf{x,y)\ < Cp{PT\frYH^,y)(^-%+ =Y T>0,{x,y)e 



Consequently, 



(1.7) ri(PT/)<^r^, T>0,/G=^, 

holds for some constant C > 0, where Ti is given by fll.3p . 

Let Pt{z; ■) be the transition probability kernel of Pt. It is easy to see that (11.60 implies 

\\Pt{{x, y)-.)- Pt{{x\ 2/0; ■) II™. < C||/|U^^ + T > 0, (x, y), {x'y') e R^^' 

for some constant C > 0, where ||v?||i,ar := sup (/?(■) — inf </)(■) is the total variational norm 
of a signed measure 93. Consequently (cf. [15]), the Markov process has successful cou- 
plings. Moreover, according to the following result, (II. 6p and (I1.7P also imply Harnack type 
inequalities for Pt- 

In general, let i? be a connected differential manifold and let Fi be a square field of type 



ri(/) = Y.(^^ff 



i=l 



for some continuous vector fields For any vector v G T^E, the intrinsic norm of v 

induced by Fi is 

l^lr, =sup{|^/|(a;):Fi(/)(x)<l}. 
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For any C^-curve 7 : [0,1] E, the length of 7 induced by Fi is 

^(7) = / l7s|rids. 
Jo 

Finally, for any z, z' G E, the intrinsic distance between them induced by Fi is 

p{z, z') = inf {£(7) : 7 is a C^-curve linking z and z']. 

It is well known that p is finite if {Xi}f^^ are smooth vector fields satisfying Hormander's 
condition. An alternative way to define p is to use the subunit curve. Recall that a C^-curve 
7 : [0,T] ^ ^ is called subunit w.r.t. Fi if \j-Jht)\ < V'^i{f)ht), t e [0,T]. Then 

p{z, z') = inf |T > : there exists a subunit curve 7 : [0, T] — M, 70 = z, = 2;'}. 

Proposition 1.3. Let Fi and p be fixed as above on a connected differential manifold E 
such that p is finite. Let P be a {sub-) Markov operator on ^fy{E), the set of all bounded 
measurable functions on E. Then for any constant C > 0, 

(1.8) T,iPf)<C'Pf, feCliE) 
is equivalent to the Harnack type inequality 

(1.9) Pfiz')<Pfiz) + Cp{z,z')y^PP{^, z,z' eEJ>OJe.^kiE). 

A simple application of (11. 9p is the following Harnack inequality for the transition kernel 
P{z, ■) of P: taking / = 1^ in (11. 9p for measurable set A, we obtain 

P{z, ■) < P{z', ■) + Cp{z, z') v/p(iV), z, z' e E. 

We will prove Theorem 11.11 in Section 2 and prove Corollary 11.21 and Proposition 11.31 in 
Section 3. Finally, in section 4 we extend Theorem II. II to a more general model. 



2 Proof of Theorem 1.1 



To establish the derivative formula, we first briefly recall the integration by parts formula 
for the Brownian motion. Let T > be fixed and let 



H 



/iG C([0,T];M™+'') : /i(0) = 0, := ^ |/i'(t)|Mt < ooj 



be the Cameron-Martin space. Let p be the distribution of {Bt, Bt)ti^[o,T], which is a prob- 
ability measure (i.e. Wiener measure) on the path space W = C([0, T]; M™+'^). A function 
F e L'^{W] p) is called differentiable if for any /i G H, the directional derivative 

e^O e 
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exists in L'^{W; /i). We write F G ^{D) if moreover 

m3h^ DhF e L'^iW] fi) 

is a bounded linear operator. In this case the Malhavin gradient DF is defined as the unique 
element in L'^(W H; /i) such that {DF,h)a = DhF for h e M. It is well known that 
{D,^{D)) is a closed operator in L^(VF;//), whose adjoint operator (5,^(5)) is called the 
divergence operator. That is, 

(2.1) / DhFdii= I F5{h)dn, F e ^iD),he ^(6). 



Theorem 2.1. For fixed T > and v = {vi,V2) G R'^+'^, let hi e C\[0,T];W^) with 
hi{0) = and hi{T) = Vi. If there exists a process {/i2(0}te[o,r] on M*^ such that /i2(0) = 0, 
and h := [hi, /i2) G ^(5) satisfying 



Proof. From (11.41) it is easy to see that the derivative process (Vt,Xt, V,,yi)t>o solve the 




(2.2) 




then 



V.Pt/ = E{/(Xt,Ft)W}, feCliR'). 



equation 




So, 



(2.3) 




Next, for h given in the theorem, we have 




dDhYt = a{Xt)h'2{t)dt + {VD,x,a){Xt)dB, 



DhXo = 
DhYo = 0. 



Thus, 



(DnXt = hi{t), 

\DhYt = J^aiXs)h'2is)ds + /o*(V,,(,)a)(X,)d5, 



Since hi(T) 



Vi, combining this with ( 12. 2p and ( 12. 3p we obtain 



(V.Xt, V.Ft) = (A.XT,I)h>T). 
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Therefore, for any / G C^(M^), it follows from ( 12. ip that 

V.Pt/ = E{Vf{XT,YT), (V.Xt, V^Ft)) = E(V/(Xt,1t), (D/.Xt, D/.Ft)) 
= ED4/(X7^, Yt)} = E{/(Xt, yT)5(/i)}. 



□ 



To prove Theorem 11.11 the key point is to solve the control problem (12. 2p . To this end, 
we will need the following fundamental lemma. 

Lemma 2.2. Let pt be a predictable process on M.'^ with E Jq \pt\'^ < oo for some q > 2. 
Then 



E 



Proof. It suffices to prove the first inequality since the second follows immediately from 
Jensen's inequality. Let Nt = jQ{ps,dBs), t>0. Then d{N)t = \pt\'^dt and 

dNf = 2NtdNt + \pt\''dt. 

Noting that |iVt|« = {N"})"/'^, by Ito's formula we obtain 

Therefore, 

EliVrl" = ^^^^ f ^{\Ntr^\pt\'']dt 
2 Jo 

<^^^ r(E|iv,r)(-^)/^(E|p,r)^/^dt 

^ Jo 

2 io 

Up to an approximation argument we may assume that E|A'r|'' < oo, so that this implies 

fT \ q/2 



E|iv.r<{^r(/ (EiP.rr/^dt)' 



□ 
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Proof of Theorem \l.l[ We assume that {Xq,Yq) = {x,y) and simply denote E^'^ by E. Let 



(2.4) 
and 

(2.5) /i2(t) 



hi{t) = ^-^, te[0,T] 



a(X,)*ds Q-M 172 + 



T-s 
T 



(V,„a)(X,)d5j, te[0,T]. 



Then it is easy to see that (12. 2p holds. To see that h := {hi, G ^{5) and to calculate 
5{h), let 

T-s ~\\ 
-(V.,a)(X,)dSj ), 



T 



9i = (^ei,Qj^^\^V2 

hiit) = / a{Xs)*eids, i = l,---,d, 
Jo 

where {ejjf^^ is the canonical ONB on M'^. We have 

d 

(2.6) h{t) = {h{t),0) + J29^i0Mt))■ 

i 

It is easy to see that hi and hi are adapted and 

5((/ii,0)) = 



i=l 



{hi{t),dBt) = — - — , 



{hr{t)ABt)= / (a(Xi)*e„dS,). 



(2.7) 



Let ^ be the cr-field induced by {Bg : s G [0,T]}. By Lemma 12.21 and noting that Xt is 
measurable w.r.t. ^ while B is independent of we have 



<2\\v2r\\Q:r'r^ 



T 



{aiXtYe^dBt) 



+ 2E 



<c\\Q^'r 



T 



T-t 



T 



{V.,a{Xt)yiQ^'re„dBt)\ { / (a(Xi)*e„ d^^) 



T 



+ 2 



E 



||a(XO||Mt 



T 



{W.,MXt)yiQT're^,dBt) 



^ -EN / (a(Xi)*e„dS4) 



1/2 



<c'iiQ?Nr 



T 



(iiva(x,)r+iia(x,)r+i)dt 



for some constants c, c' > 0. So, (A) implies giS{{0,hi)) G L^(P) for i = 1, ■ ■ ■ d. Hence, if 
for any ie {l,---d} one has ^(o,h,)5'i e ^^(P), then h G ^{6) and by ([2J]) and (ETD, 

(2.8) 5(/i) = + |j |^;.^^aM*e„dS,) - D(o,^^)^.}. 

Noting that = x + i?i is independent of it is easy to see that 

D^oM)g^ = {e,, Qt' ^ ^(V.,a)(Xi)/i:(t)dt 

= (e.,QT' ^^^{(V.,a)a*}TOe,dt 

which is in L^(P) according to (A). Combining this with (12.8P and noting that Xf = x + Bt, 
we conclude that h G ^(5) and (5(/i) = M^. Then the proof is finished by Theorem 12.11 □ 



3 Proofs of Corollary 11.21 and Proposition 11.3 



To verify (A) for cr given in Corollary II. 2[ we first present the following lemma. 
Lemma 3.1. For any n G [1, oo) and a > 0, there exists a constant c > such that 

rT 



x,y 



|XJ2"dt 1 < 



r"(|x|2 + T)^ 



-, T>0,{x,y)e 



Proof. We shall simply denote E^'^ by E. Since Xt = x + Bf, for any A > we have (see e., 
[101 page 142]) 



exp[-^tanh(V2AT)] 



.^^ {coth (V2AT) l'"/^ 



< 2^/2 exp 

< 2'"/2 exp 

< 2'"/2exp 



V2 2V2 
ix' + T)V\ 

+ T)VX ' 
2V2 



{V2XT) A l| 
{V2XT) A l| 



[X 



I 2™/2 



exp 



(x^ + T)AT 

7! . 



This implies that for any r > 0, 
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Eexp 



< Eexp 

L 

< 2™/2 exp 



n — 1 



= Eexp 









x2 + r)Ai/(2-)yf 



j^n/(n-l) 

Taking r = we obtain 



exp 



E exp 



-A TlXtpMi 
Jo 



< C\ exp 



2V2 

[x" + r)A^/(^") 



+ exp 



exp 



_ (x^ + T)rA^/" 

(x^ + T)(Ar)V" 
V2 



for some constant Ci > 0. Noting that 







e s ds - 



holds for all / > 1 and Q^a > 0, we conclude that 

A""^Eexp 



E^ |Xip"dt 



<ciy A""^<; exp 



< 



(|xp + r)A^/(^") 

2^y(n-l)/2n 
C3 



A r\xtf^dt 

Jo 



dX 



< 



+ exp 

c 



(|xp + r)(Ar)^/"' 
7! , 



dA 



holds for some constants C2, C3 and c. 



□ 



Proof of Corollary By Jensen's inequality, it suffices to prove for j» G (1,2] so that 
q := > 2. In fact, once (11. 6p holds for p = 2, it also holds for p > 2 with Cp = C2 since 

in this case {Prf^Y'^ < (PtI/I^)^/^- 
It is easy to see that (11.51) implies 



Qt> [a^ \Xtfdt^h>,d. 



and hence, 
(3.1) 



Since {Xt}teio,T] is measurable w.r.t. ^ and due to (11.51) 



2| V |2Z-1 



||{(V.,aK}(X,)||<6lX,| 
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we obtain 



E 



(3.2) 



T 

{vi,Bt) 



-1 

T 



T 



Ti[Q 



T-t 



T 
T-t 
T 



{iW,,a)a*}iXt)dt 
{i\/,,a)a*}iXt)dt 



for some constant ci > 0. Moreover, since Bt is independent of due to (13. ip and Lemma 
12.21 there exist constants C2,C3 > such that 



E( IIQ^^II'^ 



\ cj{Xt)dBt 



and 



E( iig^^r 



T-t 
T 



{V,,a){Xt)dBt, / a(Xi)d5i) 



< 



C2 



X J E 



E 



T 



T-t 



a{Xt)dBt 



T 

{2l^l)q)/l 



{y.,a){Xt)dBt 

1/(21-1) 



(2Z-l)g/a-l) 



(Z-l)/(2/-l) 



< 



hold. Combining these with (13. 2p we obtain 

E|Mt|^ = E{E(|MTn^)} 



(3.3) 



< C3E 



T^/^ ' ilolX^dty ' (/o'lX.P'dt)'^ /■ 

By Lemma [3.11 and noting that Xf = x + Bf, we conclude that for any /? > 1, 



E 



ilolX^dt)^ 



< <^E 



\Xt\^dt 



T 



< 



< 



CiT{\x 


2 _^ 2^^/3/2 


T<i{\x 


2 + 



E( / iXtl^'dt 
C4 



2g>. 1/2 



holds for some constants €3,04 > 0. Substituting this into (13.30 we arrive at 





^^2 




r'?/2(|a; 


2 _^ 2^)gV2 



1/q 
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for some constant C5 > 0. Therefore, (II .Gp follows since according to Theorem ll.il 



□ 



Proof of Proposition \1.3[ ( ll.Sp =^ ( II. Oft . By the monotone class theorem, it suffices to 
prove (11. 9p for / G Cb{E). For z, z' G i?, let p = p{z, z'). Up to an approximation argument 
we assume that p is reached by a subunit curve 7 : [0, p] — j- with 70 = 2, 7p = 2:'. Then, 
due to (11. Sp . for any positive / G Cb{E) we have 



ds V 1 + rs/ 



C2 



< 

4r 



Integrating over [0,p] w.r.t. ds we obtain 
Combining this with the fact that 



/ r rpf 



2 



/ — > / - rpf , 

1 + rpf 1 + rpf 

we obtain 

Minimizing the right-hand side in r > we prove (II. 9p . 
( ICT ^ (HSI). By dini), we have 

|P/(x)-P/(zOI <Cp(x,?/)||/|U, /ea(M). 

So, P/ is p-Lipschitz continuous for any / G ^b{E). Let z E E and 7 : [0, 1] — )■ M be 
C^-curve such that ■jo = z, p(7o, 7s) = s and 

;^P/(7.)|.=o= v/r(P/)(.). 
ds 

Then it follows from (II. 9p that 

yf(P7M = lim ^^^^'^ ~ < C hm v/PfM = C^/Ppiz). 

Therefore, ([LlD holds. □ 
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4 An extension 

Consider the following SDE on ]R™'+'^ 
(4.1) 



dXt = ai{Xt)dBt + h{Xt)dt, 
dYt = a2{Xt)dBt + b2{Xt)dt, 



where {Bt, Bt) is a Brownian motion on ai G Cfe^(M";]R™ M™) is invertible with 

||o"r^ll — some constant c > 0, G C^(M'";M'^ (g> IR'^) might be degenerate, bi G 

Cfe^(M™; M") and 62 G C^(M™; M'^). It is easy to see that for any initial data the solution exists 
uniquely and is non-explosive. Let Pt be the associated Markov semigroup. To establish the 
derivative formula, let v = {vi, V2) G MJ^~^'^ and T > be fixed, and let solve the following 
SDE on M'": 

(4.2) det = (V5,ai)(X,)di?, + {(V5,MW-^}dt, 6 = 0. 

Since Va and V61 are bounded, the equation has a unique solution up to time T. It is easy 
to see from the Ito formula that 



T-tj V \^ T-t T-t (T-t) 

< 2(^, {V^,a){XMBt) + - (^ri^y^jd^, t G [0,T) 

holds for some constant C > 0. This implies that for t G [0, T), 

(4.3) E|e<P<(T-t)e^*, E^'^J^dt<oo, 

Consequently, we may set = so that solves (14. 2 p for t G [0,T]. Moreover, for any 
n > 1 we have 

dietl'" < 2n|6r("-^)(6, iV^,a){Xt)dBt) + c{nm'''dt 
for some constant c(n) > 0. Therefore, 

(4.4) sup EI^T" < 00, n > 1. 

ie[o,T] 



We are now able to state the derivative formula for Pf as follows. 
Theorem 4.1. Let Qt = (72{Xt)a2{Xt)*dt he invertible such that 



(4.5) 



W'y[\\Q^^r j\\\a2{Xt)\\' + ||Va2(X0r + \\Vh2{Xt)f + l}dt) < 00. 



13 



Then 

V,PTf{x,y) = E-'y{fiXT,YT)MT} 
holds for f G Cl{W^+'^) and 



Mt=1 (^^^r^'di?.)-Tr(g^^^ ^{(V,,a,)a;}TOdt) 

Proof. Let h = {hi, /i2), where 
Jo T - s 

As in the proof of Theorem ll.il it is easy to see from gSD, diaD, (531) and Wa^'^W < cthat h G 
^(5) with = My. Therefore, it remains to verify that (V^Xt, V^Yt) = {DhXr, DhYr). 
It is easy to see that both V^Xf and D^Xf + C,t solve the equation 

dVt = {Vv,ai){Xt)dBt + {VvM)iXt)dt, t G [0,T],Vo = v^. 

By the uniqueness of the solution we have VyXt = DhXt + for t G [0,T]. Since = 0, 
this implies that V^Xy = DhXr- Moreover, we have 

'dV.Fi = {V^,,x,a2){Xt)dBt + {V^^xM){Xt)dt, V.Fo = V2, 
dDhYt = {VD,x,a2)iXt)dBt + a2iXt)h'^{t)dt + {V D,xM)iXt)dt, DhY^ = 0. 

Combining this with the definition of /i2 and D^Xt = V^X^ — ^t, we obtain 

DhYT = V,YT-V2- [ {V^,a2){Xt)dBt+ [ a2{Xt)h'^{t)dt - f {V^M){Xt)dt = V,Yt. 

Jo Jo Jo 

Therefore, the proof is finished. □ 

Acknowledgement. The author would like to thank the referee for careful reading and 
useful comments. 
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